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In this work, we suggest that hard-pion chiral perturbation theory may be appHcable to the real 
parts of nonleptonic — > D^P^ and T^P^ (P = k) decay amplitudes. These amplitudes 
play an important role in the extraction of the angle y in the b—d unitarity triangle of the CKM 
matrix, and their real parts can be computed using lattice QCD. We construct the leading-order 
operator in the chiral expansion for these nonleptonic decays, and discuss the generic features of 
the next-to-leading-order terms. 
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1. Introduction 

Nonleptonic B decays have played an important role in CKM physics. In particular, the angle y 
in the b—d unitarity triangle can in principle be determined with high precision from the study of 
charged B meson decays processes, 

B+ ^ D*'P+ ^ /P+ and B+ D°P+ ^ /P+, (1.1) 

with the same final state /P+ [|l|, ^ ||, ^]. In all the methods for extracting 7 from the above 
decay channels, the knowledge of the branching ratios, Br[B^ — )■ D°P+] and Br[B^ — D^P^], is 
a key ingredient. While Br[B^ — D^P+J has been experimentally measured with good accuracy, 
Br[B+ is very difficult to obtain. For this reason, despite the large statistics of the two 

B-factories [~ ^(10^) charged B meson samples], 7 is presently determined to only ~ ^(25%). 
This is compared with about 3% for j3, and about 5% for a. To further improve the precision on 7, 
inputs from lattice QCD (LQCD) to these branching ratios will be very helpful 

The study of hadronic weak decays on the lattice is very challenging, because of the Maiani-Testa 
no-go theorem (MTNGT) [Q]. This theorem states that by using Euclidean four-point correla- 
tors to study nonleptonic two-body decays, it is impossible to obtain information about the strong 
phases. Therefore one can only compute the real parts of nonleptonic decay amplitudes from such 
correlators. For the calculation of K ^ nn on the lattice, one can avoid the MTNGT using the 
Lellouch-Liischer (LL) method [^. On the other hand, the lattice computation of nonleptonic B 
decays remains challenging, because the LL method is only applicable to processes involving elas- 
tic final-state scatterings. Nevertheless, lattice results for the real part of these amplitudes could 
provide valuable information of Br[B^ — )• D°P+] and Br[S+ — )• D*'p+], and help in the extraction 
of 7. 

In the near future, lattice computations for such four-point functions are unlikely to be preformed 
exactly at physical pion mass. Thus, it is important to rely on the chiral expansion to perform 



extrapolations to the physical point. In Ref. 1 15], we examine the possibility of the chiral expansion 
of the real parts of B+ — > D°P+ and B+ — >^ D P+ amplitudes, in the framework of heavy-meson 



chiral perturbation theory (HMZPT). [[7|, |8|, |9|, |10[ |ll|]. The straightforward application of the chiral 



expansion for B — DP amplitudes near the physical kinematics is questionable, since the final-state 
hadrons carry momenta ~ 2 GeV. From the structure of subleading terms in l/M^ (Mp is the D 
meson mass) in the effective theory [12, 13, |l4|, it is clear that the large momentum carried by the D 



meson will only lead to (significant) dependence on Md in the low-energy constants (LLC's). What 
remains to be considered is the convergence of the chiral expansion for hard final-state hadrons 
from lattice studies. 



In Ref. [|15|], we address this issue in nonleptonic B decays, in a framework that is similar to the 
hard-pion chiral perturbation theory (HP;^PT) which was first proposed and applied to the analysis 
of Kgj, decays in Ref. [|l6|]. This approach was also used to carried out investigations of K ^2n 
amplitudes [p^^ , as well as semi-leptonic B-decays [18]. 



' As pointed out in Ref. hard-pion rnay not be used to estimate the strong phases in A" — > JlK processes. 
We will comment on this issue for nonleptonic B decays in Sec. H. 
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The relevant current-current, AZ? = 1, operators (a,j3 are colour indices) in our work are, 

e^''' = (^L7^1-75)^«)(c/37m(1- 75)^/3) , (1-2) 
e'^''' = (^L7^(l - 75)bii){cpy^{l - 75)"a) , (1.3) 
e^"''' = (^^7^(1 - 75)^«)(«/37m(1 - 75)c/3) , (1.4) 

Ql"""' = (^a7^(l - 75)^/5)(%7m(1 - 75)Ca) • (1.5) 

We will focus on the nonleptonic decays with the underlying processes b — cud, b — cus, b — ucd, 
and b — )■ mc5. The first two will be mapped onto different operators in the chiral effective theory 
from the last two, as explained in the next section. 

2. The leading-order chiral expansion for B — DK{n) amplitudes 

In weak nonleptonic B decays, the final-state hadrons carry large momenta. This makes it difficult 
to apply the chiral expansion to such processes. On the other hand, it has been recently proposed 
that chiral perturbation theory (;tPT) can be valid for amplitudes containing hard final-state par- 
ticles [|^, |l^, jTsj, |l9|]. One crucial point is to allow the low-energy constants (LECs) to depend 
on the hard momentum scales which result from either the kinematics or the mass of the external 
particles. That is, the LECs in the chiral expansion are no longer universal quantities. Another 
important ingredient in such chiral expansions is the separation of the hard and soft scales. This 
separation of scales is made possible because of the structure of derivative couplings in ;tPT, lead- 
ing to the absorption of the hard, external, momenta into the LECs. We will discuss this procedure 
explicitly for B — )• DP amplitudes with an example diagram in Sec. ^ 



In this section, we construct the XVT weak operators corresponding to those in Eqs. (|1.2|)-(|1.5|) 



Neglecting the colour indices which do not play a role in ZPT, these operators can be written as 

e^^^'' = (^iri b) (cr2«L), 



where q' =d ox s, and 



Q'^'^' = {qir^b){uLr2c), (2.1) 



qL={Kr^]q; ri=r2 = ri=r2 = 7M(i-75)- (2.2) 



Under the SU(3)l (E> SU(3)r chiral symmetry group, g^^^ ' is in the (8l, 1r) representation, while 
Qb^c,i jj^g (6l,1r) representation. To bosonise these operators, we promote Fi 2 and Fi 2 to 
be spurion fields which transform as 

Fl^LFlS^ F2^5F2L^ ^l^Lfl5^ ^2^Lf2 5^ (2.3) 



under the heavy-quark spin/flavour and chiral rotations. This renders the operators in Eq. (2.1) 



invariant with respect to such transformations. Using this property, and taking into account all the 
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possible insertions of Dirac structures [pO|], we obtain the leading-order (LO) operators, 



X,' 



i8i+i82(v'-v) 



4^2-^3-2 i8i+i84 + i85 (v'.v 



(2.4) 



where j8,- and j8,- are LECs, v and v' are the velocities of B and D mesons, a = a/T =exp(/<I>//) 
with <I> and £ being the standard linear and nonlinear Goldstone fields, ^C*) and y*^'^^ are the 
anihilation fields for the pseudoscalar and vector mesons containing the heavy quark h. At the LO, 



(2.5) 



This is not valid at the next-to-leading-order (NLO) chiral expansion. 
3. Resonance contributions 

In this section we investigate resonance contribution to B — )• DP correlators (in the time-momentum 
representation) and amplitudes^. The vector heavy-light resonances are already incorporated in 
HM;fPT. The inclusion of heavier resonances in the effective theory is beyond the scope of this 
work. Because of the complications in formulating HMZPT in Euclidean space, and here we are 
only sutdying the generic feature of the resonance contribution, we work in Minkowski space. 



^The conclusion presented in this section is also valid for B DP decays. 
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B ' D B ' ^ F 

D' 

(a) (b) 

Figure 1: Tree-level correlators contributing to B DP. The box is the weak operator, with (a) being the 
direct B DP term and (b) being the term with an intermediate resonance (here a D*). 

First, we calculate the LO correlator in Fig. [l](a). We perform a Fourier transform for the spatial 
directions for each of the external-source points, while fixing the location of the weak operator (the 
square in the diagram) to be at the origin. As in Sec. ^, we denote the velocity of B and B* by v and 
that of D and D* by v'. For simplicity, we choose the frame in which v = implement 
the time-ordering as Ib < < to < tp, where tB.D.p are the temporal locations of the sources for B, 



D and P mesons. Using the the weak operators in Eq. (2.4), the result of this diagram is 



Clo 




(3.1) 



where 5d = v' ■ po, and (Op = \ Mp + pp, with po and pp denoting the spatial momenta of the D 



and the P mesons. The coupling gBUP is one of the linear combinations of the LEC's j8, in Eq. (2.4) 
Next, we compute the diagram in Fig. [T](b), which leads to the result 



gBD*gn 




Qi{(Op + SD-ADp)tD _ I 

2v'{(Op + 5d-^dp) 



(3.2) 



where Ap,p = v' ■ {po + Pk) + Ad/v'q , with Ao denoting the D*—D mass splitting resulting from the 
heavy-quark spin symmetry breaking effects. The coupling gji is the LO B*—B—n axial coupling in 



the HM;^PT Lagrangian, and gBo* is a linear combination of the LEC's j8, in Eq. {2A). Notice that 
gBD' is different from gBOP and thus the resonance contribution results in general in an additional 
unknown parameter for B — )• DP amplitude at the tree level. When the final-state momenta are tuned 



such that the resonance is on-shell, the factor in the square brackets in Eq. ( |3.2| ) becomes linear in 
to- This can be interpreted as an energy shift of the final state. In any case, the to dependence in 
these square brackets will be cancelled when taking the ration between Clo + Qes and the square 
root of the DP — )• DP correlator to obtain the B — )• DK matrix element. 



4. Hard-pion chiral perturbation theory at one-loop 

To account for nonanalytic dependence on niyc in B — )■ DK amplitudes, it is necessary to go beyond 
the tree-level in the chiral expansion. In order to treat these processes in the physical regime, we 
use the methods of Refs. [p^, 17, 18]: Hard-pion ZPT (HPZPT), in which one key point is the 



separation of the hard external momenta with the soft scales (the Goldstone masses) which appear 
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Figure 2: One of the many one-loop diagrams that contribute to B — > Dn, specifically one which shows the 
essential features that arise in HPZPT. 

in the loop. This is acheivable because of the structure of the derivative couplings. Focusing on the 
SU(2) xVT, the chiral logarithms for B — )• DK amplitudes in HP^PT takes the generic form, 



2 / 2 



(4.1) 



where ^'^'^^'^ is its LO value, a is a coefficient that depends on the particular kinematics chosen for 
the diagram, and L is a linear combination of low-energy constants as well as terms arising from 
higher-order chiral-level weak operators. The coefficient a would be determined from evaluating 
the one-loop corrections to the amplitude. Both a and L depend on all of the hard quantities (such 
as the momenta of the external D meson and pion). 

In order to understand the specific details, we examine the diagram in Fig. ^, through which we 
will demonstrate the dependence on the hard scales in the coefficient of the chiral logarithm. In this 
diagram, as in lattice calculations, momentum will be conserved at the strong vertex, but need not 
be at the weak operator. This momentum insertion, p^\^, is related to the momenta of the external 
mesons via + Pwk = 7^0 + P%- Denoting the residual momentum of the external D by k, and the 
internal pion momentum by /, the integral that needs to be computed for this diagram is 



dH 



IV 



■Pn) 



(271)4 ^2_„^2+j-£^/.(^_^_^^)_^ + j-£ ' 



(4.2) 



where v' is the velocity of the D mesons, and A = mi, — tub is the D-B meson mass splitting. This 
integral can be evaluated using dimensional regularisation. 



v'-k + A 



v'.(^ + ;..)+A + /s-'^("-^'-('^ + ^^) + ^ + '^)-'"'^l'^(S 



(4.3) 



where h takes the form in the hard-pion limit v'-k » irijc, 



l2{m^,v'ik + p^) + A)^-ml\n{-^] , 



mt 



(4.4) 



so that the full coefficient a [as defined in Eq. (l-.l) ] in the integral / is —2 when we insert mo- 
mentum into the weak vertex such that pj^ 0, and is —3/2 when we choose ;?wk such that pji k. 
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Here we comment that the imaginary part in this diagram is proportional to \/ [v' • {pji +k)Y — m\, 
therefore grows with the final-state momenta, leading to the failure of the chiral expansion when 
Pb and pji are large. This can be understood by noting that the imaginary part arises from the 
contribution in which both mesons in the loop are on-shell, and therefore cannot be soft. 

5. Summary 

We studied the chiral expansion for lattice computations of the real parts of the — D^P^ and 
S+ — )• D*'/'+, which are important inputs in the determination of the angle 7 in the b—d untarity 
triangle in the CKM matrix. The calculation of these real parts are not obstacled by the MTNGT. 
We derived the LO operators for these matrix elements in HM;|^PT, argue that HP;|^PT is appUcable 
to these computations, and investigate some generic features of the NLO contributions. 
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